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Abstract 

The dynamics of many social, technological and economic phenomena are driven by 
individual human actions, turning the quantitative understanding of human behavior 
into a central question of modern science. Current models of human dynamics, used 
from risk assessment to communications, assume that human actions are randomly 
distributed in time and thus well approximated by Poisson processes . In con- 

trast, there is increasing evidence that the timing of many human activities, ranging 
from communication to entertainment and work patterns, follow non-Poisson statis- 
tics, characterized by bursts of rapidly occurring events separated by long periods of 
inactivity |^[ 5, (], 7, Here we show that the burst y n ature of human behavior is 
a consequence of a decision based queuing process |i3. llo|: when individuals execute 
tasks based on some perceived priority, the timing of the tasks will be heavy tailed, 
most tasks being rapidly executed, while a few experience very long waiting times. In 
contrast, priority blind execution is well approximated by uniform interevent statis- 
tics. These findings have important implications from resource management to service 
allocation in both communications and retail. 
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Humans participate on a daily basis in a large number of distinct activities, ranging 
from electronic communication, such as sending emails or making phone calls, to browsing 
the web, initiating financial transactions, or engaging in entertainment and sports. Given 
the number of factors that determine the timing of each action, ranging from work and 
sleep patterns to resource availability, it appears impossible to seek regularities in human 
dynamics, apart from the obvious daily and seasonal periodicities. Therefore, in contrast 
with the accurate predictive tools common in physical sciences, forecasting human and social 
patterns remains a difficult and often elusive goal. 

Current models of human activity are based on Poisson processes, and assume that in 
a dt time interval an individual (agent) engages in a specific action with probability qdt, 
where q is the overall frequency of the monitored activity. This model predicts that the 
time interval between two consecutive actions by the same individual, called the waiting or 
inter-event time, follows an exponential distribution (Fig. ^ a-c) fl]. Poisson processes are 
widely used to quantify the consequences of human actions, such as modelling traffic flow 
patterns or accident frequencies [l], and are commercially used in call center staffing Q], 
inventory control |3j, or to estimate the number of congestion caused blocked calls in mobile 
communications 4]. Yet, an increasing number of recent measurements indicate that the 
timing of many human actions systematically deviate from the Poisson prediction, the wait- 
ing or inter-event times being better approximated by a heavy tailed or Pareto distribution 
(Fig. ^ d-f). The differences between Poisson and heavy tailed behavior is striking: a Pois- 
son distribution decreases exponentially, forcing the consecutive events to follow each other 
at relatively regular time intervals and forbidding very long waiting times. In contrast, the 
slowly decaying heavy tailed processes allow for very long periods of inactivity that separate 
bursts of intensive activity (Fig. 1). 

To provide direct evidence for non-Poisson activity patterns in individual human behav- 



ior, we study the communication between several thousand emai 



users based on a dataset 



capturing the sender, recipient, time and size of each email |lll.ll2|. As Figure 2a shows, the 
distribution of the time differences between consecutive emails sent by a selected user is best 
approximated with P(r) ~ t~°, where a ~ 1, indicating that an individual's email pattern 
has a bursty non-Poisson character: during a single session a user sends out several emails in 
a quick succession, followed by long periods of no email activity. This behavior is not limited 
to email communications. Measurements capturing the distribution of the time differences 
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between consecutive instant messages sent by individuals during online chats [5] show a sim- 
ilar pattern. Professional tasks, such as the timing of job submissions on a supercomputer 
[f|, directory listings and file transfers (FTP requests') initiated by individual users or 
the timing of printing jobs submitted by users |13] were also reported to display non-Poisson 
features. Similar patterns emerge in the time interval distribution between individual trades 
in currency futures jq| . Finally, heavy tailed distributions characterize entertainment related 
events, such as the time intervals between consecutive online games played by the same user 

0. 

The fact that a wide range of human activity patterns follow non-Poisson statistics sug- 
gests that the observed bursty character reflects some fundamental and potentially generic 
feature of human dynamics. Yet, the mechanism responsible for these striking non-random 
features remain unknown. Here we show that the bursty nature of human dynamics is a 
consequence of a queuing process driven by human decision making: whenever an individual 
is presented with multiple tasks and chooses among them based on some perceived priority 
parameter, the waiting time of the various tasks will be Pareto distributed. In contrast, first- 
come-first-serve and random task execution, common in most service oriented or computer 
driven environments, lead to a uniform Poisson like dynamics. 

Most human initiated events require an individual to weigh and prioritize different activ- 
ities. For example, at the end of each activity an individual needs to decide what to do next: 
send an email, do some shopping, or place a phone call, allocating time and resources for 
the chosen activity. Consider an agent operating with a priority list of L tasks. After a task 
is executed, it is removed from the list, offering the opportunity to add another task. The 
agent assigns to each task a priority parameter x, which allows it to compare the urgency 
of the different tasks on the list. The question is, how long will a given task have to wait 
before it is executed. The answer depends on the method the agent uses to choose the task 
to be executed next. In this respect three selection protocols are particularly relevant 
for human dynamics: 

(i) The simplest is the first-in-first-out protocol, executing the tasks in the order they 
were added to the list. This is common in service oriented processes, like the first-come- 
first-serve execution of orders in a restaurant or getting help from directory assistance and 
consumer support. The time period an item stays on the list before execution is determined 
by the cumulative time required to perform all tasks added to the list before it. If the time 
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necessary to perform the individual tasks are chosen from a bounded distribution (i.e. the 
second moment of the distribution is finite), then the waiting time distribution will develop 
an exponential tail, indicating that most tasks experience approximately the same waiting 
time. 

(ii) The second possibility is to execute the tasks in a random order, irrespective of their 
priority or time spent on the list. This is common, for example, in educational settings, when 
students are called on randomly, and in some packet routing protocols. The waiting time 
distribution of the individual tasks (i.e. the time between two calls on the same student) in 
this case is also exponential. 

(iii) In most human initiated activities task selection is not random, but the individual 
executes the highest priority item on its list. The resulting execution dynamics is quite 
different from (i) and (ii): high priority tasks will be executed soon after their addition to 
the list, while low priority items will have to wait until all higher priority tasks are cleared, 
forcing them to stay on the list for considerable time intervals. In the following we show 
that this selection mechanism, practiced by humans on a daily basis, is the likely source of 
the fat tails observed in human initiated processes. 

We assume that an individual has a priority list with L tasks, each task being assigned 
a priority parameter Xi, i — 1, L, chosen from a p(x) distribution. At each time step the 
agent selects the highest priority task from the list and executes it, removing it from the 
list. At that moment a new task is added to the list, its priority X{ being again chosen from 
p(x). This simple model ignores the possibility that the agent occasionally selects a low 
priority item for execution before all higher priority items are done, common, for example, 
for tasks with deadlines. This can be incorporated by assuming that the agent executes the 
highest priority item with probability p, and with probability 1 — p executes a randomly 
selected task, independent of its priority. Thus the p — > 1 limit of the model describes the 
deterministic (iii) protocol, when always the highest priority task is chosen for execution, 
while p — > corresponds to the random choice protocol discussed in (ii). 

To establish that this priority list model can account for the observed fat tailed interevent 
time distribution, we first studied its dynamics numerically with priorities chosen from a 
uniform distribution x^ G [0,1]. Computer simulations show that in the p — > 1 limit the 
probability that a task spends r time on the list has a power law tail with exponent a = 1 
(Fig 3a), in agreement with the exponent obtained in email communications (Fig 2a). In 
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the p — > limit P(t) follows an exponential distribution (Fig. 3b), as expected for the 
case (ii). As the typical length of the priority list differs from individual to individual, it is 
particularly important for the tail of P(r) to be independent of L. Numerical simulations 
indicate that this is indeed the case: changes in L do not affect the scaling of P{t). The fact 
that the scaling holds for L = 2 indicates that it is not necessary to have a long priority list: 
as long as individuals balance at least two tasks, a bursty heavy tailed interevent dynamics 
will emerge. 

To determine the tail of P{t) analytically we consider a stochastic version of the model 
in which the probability to choose a task with priority x for execution in a unit time is 
n(x) ~ x 7 , where 7 is a parameter that allows us to interpolate between the random choice 
limit (ii) (7 = 0, p — 0) and the deterministic case, when always the highest priority item 
is chosen for execution (iii) (7 = 00, p — 1). Note that this parameterization captures 
the scaling of the model only in the p — > and p — > 1 limits, but not for intermediate p 
values, thus it is chosen only for mathematical convenience. The probability that a task 
with priority x waits a time interval t before execution is f(x,t) = (1 — U(x)) t ~ 1 U(x) . The 
average waiting time of a task with priority x is obtained by averaging over t weighted with 
f(x, t), providing 

00 11 

i.e. the higher an item's priority, the shorter is the average time it waits before execution. 
To calculate P(t) we use the fact that the priorities are chosen from the p(x) distribution, 
i.e. p(x)dx = P(r)dT, which gives 

p(r- lh ) 

In the 7 — > 00 limit, which converges to the strictly priority based deterministic choice 
(p = 1) in the model, Eq. (j2J) predicts P(t) ~ r _1 , in agreement with the numerical results 
(Fig 3a), as well as the empirical data on the email interarrival times (Fig 2a). In the 7 = 
(p = 0) limit t(x) is independent of x, thus P(t) converges to an exponential distribution, 
as shown in Fig. 3b (see Supplementary Information). 

The apparent dependence of -P(t) on the p(x) distribution from which the agent chooses 
the priorities may appear to represent a potential problem, as assigning priorities is a sub- 
jective process, each individual being characterized by its own p(x) distribution. According 
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to Eq. (0), however, in the 7 — > 00 limit P(r) is independent of Indeed, in the deter- 
ministic limit the uniform p(x) can be transformed into an arbitrary d (x) with a parameter 
change, without altering the order in which the tasks are executed llpj . This insensitivity 
of the tail to p(x) explains why, despite the diversity of human actions, encompassing both 
professional and personal priorities, most decision driven processes develop a heavy tail. 

To obtain empirical evidence for the validity of the proposed queuing mechanism we 
consider the email activity pattern of an individual Once in front of the computer, 

an individual will reply immediately to a high priority message, while placing the less urgent 
or the more difficult ones on its priority list to compete with other non-email activities. We 
propose, therefore, that the observed interevent time distribution is in fact rooted in the 
uneven waiting times experienced by different tasks. To test this hypothesis we need to 
determine directly the waiting time for each task. In the email dataset we have the time, 
sender and recipient of each email transmitted over several months by each user, thus we 
can determine the time it takes for a user to reply to a received message [ljj. As Fig. 2b 
shows, we find that the waiting time distribution P(r w ) for the user whose P(t) is shown 
in Fig. 2a is best approximated by P(t w ) ~ T~ aw with exponent a w = 1, supporting our 
hypothesis that the heavy tailed waiting time distribution drives the observed bursty email 
activity patterns. 

As in the p — > 1 limit of the model the priority list is dominated by low priority tasks, 
new tasks will often be executed immediately. This results in a peak at P(t = l)(see Fig. 
3 in the Supplementary Information), which, while in some cases may represent a model 
artifact, in the email context is not unrealistic: most emails are either deleted right away 
(which is one kind of task execution), or are immediately replied to. Only the more difficult 
or time consuming tasks will queue on the priority list. The email dataset does not allow 
us to resolve this peak, however: a message to which the user replies right away will appear 
to have some waiting time, given the delay between the arrival of the message and the time 
the user has a chance to read it. 

While we illustrated the queuing process on emails, in general the model is better suited to 
capture the competition between different kinds of activities an individual is engaged in, i.e. 
the switching between various work, entertainment and communication events. Indeed, most 
datasets displaying heavy tailed interevent times in a specific activity reflect the outcome 
of the competition between tasks of different nature. For example, the starting of an online 
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gaming session implies that all higher priority work and entertainment related activities 
have been already executed. 

Detailed models of human activity require us to consider the impact of a number of 
additional mechanisms on the queuing process. First, in the priority list model we assumed 
that the time necessary to execute a task (service time) is the same for all tasks. The size 



distribution of emails is heavy tailed 



16L Il7| , however, thus the waiting time distribution 



could be driven entirely by the time it takes to read an email, i.e. the message size. Yet, as 
Fig. 2c shows, we fail to find a correlation between the size of the email received by a user 
and the time the user takes to reply to it. While a detailed analysis should also consider 
the role of attachments, Fig. 2c suggests that the priority of a response is more important 
than the message size. Furthermore, the priorities assigned to tasks are often driven by 
optimization processes, as agents aim to maximize profits or minimize the overall time spent 
on some activity. 

A natural extension of the model is to assume that tasks arrive at a rate A and are 
executed at a rate /i, allowing the length of the priority, list L to change in time. In this case 
the model maps into Cobham's priority queue model Q, which has a power law distributed 
waiting time with a = 3/2 only when A = fi (see Supplementary Information). Thus to 
account for the power law waiting times the model requires an additional mechanism that 
guarantees A = /i (which, as Fig. 3d indicates, is not satisfied for most email users). In 
contrast, in the proposed priority list model we assumed that for humans the length of the 
priority list remains relatively unchanged (i.e. L is constant). To understand the origin 
of this assumption we must realize that for A = \i the length of the priority list fluctuates 
widely and can occasionally grow very long. While keeping track of a long priority list is 
not a problem for a computer, it is well established that the immediate memory of humans 

we can easily remember, 

and therefore the length of the priority list, is bounded, motivating our choice of a finite L. 

While other generalizations are possible and often required, our main finding is that the 
observed fat tailed activity distributions can be explained by a simple hypothesis: humans 
execute their tasks based on some perceived priority, setting up queues that generate very 
uneven waiting time distributions for different tasks. In this respect the proposed priority 
list model represents only a minimal framework that allows us to demonstrate the poten- 
tial origin of the heavy tailed activity patterns, and offers room for further extensions to 
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capture more complex human behavior. As the exponent of the tail could depend on the 
details of the prioritizing process, future work may allow the empirical data to discriminate 
between different queuing hy pot heses. A mapping into punctuated equilibrium models (see 
Supplementary Information with the mathematical framework of queuing theory 

could help the systematic classification of the various temporal patterns generated by human 
behavior. 



reavy-tailed 
As 



There is overwhelming evidence that Internet traffic is characterized by 
statistics j^, rooted in the Pareto size distribution of the transmitted files 
we have shown above (Fig 2c), we find that a user's email activity does not correlate with 
the email size. Similarly, the timing of online games ^| or sending an instant message ^ 
cannot be driven by file sizes either. This suggests that Internet traffic is in fact driven by 
two separate processes: The heavy tailed size distribution of the files sent by the users and 
the human decision driven timing of various Internet mediated activities individuals engage 
in. In some environments this second mechanism, whose origin is addressed in this paper, 
can be just as important as the much investigated first one. Given the differences in routing 
performance under Poisson and Pareto arrival time distributions 0, 0, Q> a queuing 
based model of human-driven arrival times could contribute to a better understanding of 
Internet traffic as well. 

Uncovering the mechanisms governing the timing of various human activities has signifi- 
cant scientific and commercial potential. First, models of human behavior are indispensable 
for large-scale models of social organization, ranging from detailed urban models Q E|, 
to modeling the spread of epidemics and viruses, the development of panic |2|| or captur- 
ing financial market behavior 2(|. Understanding the origin of the non- Poisson nature of 
human dynamics could fundamentally alter the dynamical conclusions these models offer. 
Second, models of human behavior are crucial for better resource allocation and pricing 
plans for phone companies, to improve inventory and service allocation in both online and 
brick-and-mortal retail, and potentially to understand the bursts of ideas and memes emerg- 
ing in communication and publication patterns j^f. Finally, heavy tails have been observed 
in the foraging patterns of birds as well 3], raising the intriguing possibility that animals 
also utilize some evolution-encoded priority based queuing mechanisms to decide between 
competing tasks, like caring for offsprings, gathering food, or fighting off predators. 



S 



Supplementary Information accompanies the paper on www.nature.com/nature 
Acknowledgements I have greatly benefited from discussions with Alexei Vazquez on 
the mathematical aspects of the model. I also thank Luis Amaral, Zoltan Dezso, Plamen 
Ivanov, Janet Kelley, Janos Kertesz, Adilson Motter, Maya Paczuski, Kim Sneppen, Tamas 
Vicsek, Ward Whitt, and Eduardo Zambrano for useful discussions and comments on the 
manuscript, Jean-Pierre Eckmann for providing the email database and Suzanne Aleva for 
assisting me with the manuscript preparation. This research was supported by NSF DMR04- 
26737 and NSF DMS04-41089. 

Competing interests statement The author declares that he has no competing finan- 
cial interests. 

Correspondence and requests for materials should be addressed to A.-L. B. 
(alb@nd.edu). 



9 



[1] Haight, F. A., Handbook of the Poisson Distribution (Wiley, New York, 1967). 
[2] Reynolds, P., Call Center Staffing (The Call Center School Press, Lebanon, TN, 2003). 
[3] Greene, J. H., Production and Inventory Control Handbook (McGraw-Hill, New York, 3 ed, 
1997). 

[4] Anderson, H.R., Fixed Broadband Wireless System Design (Wiley, New York, 2003). 

[5] Dewes, C, Wichmann, A., & Feldman, A., An Analysis of Internet Chat Systems, Proc. 2003 

ACM SIGCOMM Conf. on Internet Measurement (IMC-03), ACM Press, New York, 2003). 
[6] Kleban, S.D. & Clearwater, S.H., Hierarchical Dynamics, Interarrival Times and Performance, 

Proc. of SC'03, November 15-21, 2003, Phonenix, AZ, USA. 
[7] Paxson, V. & Floyd, S., Wide- Area Traffic: The Failure of Poisson Modeling, IEEE/ACM 

Tansactions in Networking 3, 226 (1995). 
[8] Masoliver,J. , Montero, M. & Weiss, G.H., Continuous-time random- walk model for financial 

distributions, Phys. Rev. E67, 021112 (2003). 
[9] Cobham, A., Priority Assignment in Waiting Line Problems, J. Oper. Res. Soc. Amer. 2, 

70-76 (1954). 

[10] Cohen, J. W., The Single Server Queue (North Holland, Amsterdam, 1969). 

[11] Eckmann,J-P., Moses, E. & Sergi, D., Entropy of dialogues creates coherent structures in 

e-mail traffic, Proc. Nat. Acad. Sci. USA 101, 14333-14337 (2004). 
[12] Ebel, H., Mielsch,L.-L, & Bornholdt, S., Scale-free topology of e-mail networks, Phys. Rev. E 

66, R35103 (2002). 

[13] Harder, U., & Paczuski, M., Correlated dynamics in human printing behavior, 
|http://xxx.lanl.gov/abs/cs.PF/0412027| 

[14] Henderson, T., & Nhatti, S., Modelling user behavior in networked games, Proc. ACM Mul- 
timedia 2001, Ottawa, Canada, 212-220, 30 September-5 October (2001). 

[15] Miller, G.A., The Magical Number Seven, Plus or Minus Two: Some Limits on Our Capacity 
for Processing Information, Psychological Review, 63, 8197 (1956). 

[16] Crovella, M. & Bestavros, A., Self-similarity in World Wide Web traffic: evidence and possible 
causes, IEEE/ACM Trans. Netw. 5, 835-846 (1997). 

[17] Mitzenmacher, M., A brief history of generative models for power law and lognormal distri- 



10 



butions, Internet Mathematics 1, 226-251 (2004). 
[18] Bak, P. & Sneppen, K., Punctuated equilibrium and criticality in a simple model of evolution, 

Phys. Rev. Lett. 71 4083-4086 (1993). 
[19] Jensen, H.J. Self- Organised Criticality (Cambridge Univerity Press, Cambridge, 1998). 
[20] Park, K. & Willinger, W. Self-similar networks traffic and performance evaluation (Wiley, 

New York, 2000). 

[21] Leighton, F.T., Maggs, B.M., & Rao, S.B., Packet routing and job-shop scheduling in 

0(congestion + dilation) steps, Combinatorica, 14 167-186 (1994). 
[22] Harris, CM., Brill, P.H., Fischer, M.J. Internet-Type Queues with Power- Tailed Interarrival 

Times and Computational Methods for Their Analysis, INFORMS J. on Computing 12, 261- 

271 (2000). 

[23] Eubank, H., Guclu, S., Kumar, V.S.A., Marathe, M., Srinivasan, A., Toroczkai, Z., & Wang, 
N., Controlling Epidemics in Realistic Urban Social Networks, Nature 429, 180-184 (2004). 

[24] Manrubia,S. C, Zanette, D., & Sole, R. V., Transient dynamics and scaling phenomena in 
urban growth, Fractals 7, 1-8 (1999). 

[25] Helbing, D., Farkas, I. & Vicsek, T., Simulating Dynamical Features of Escape Panic, Nature 
407, 487-490 (2000). 

[26] Caldarelli, C, Marsili, M. & Zhang, Y.-C, A Prototype Model of Stock Exchange, Europhysics 

Lett. 40, 479-484 (1997). 
[27] Kleinberg, J., Bursty and Hierarchical Structure in Streams, Proc. of the 8th ACM SIGKDD 

Int'l. Conf. on Knowledge Discovery and Data Mining, 91-101 (2002). 
[28] Viswanathan, G. M., Buldyrev, S. V., Havlin, S., da Luz,M. G. E. , Raposo,E. P., & Stanley, 

H.E., Optimizing the success of random searches, Nature 401, 911-914 (1999). 



11 



500 





200 400 600 

Event Number 



20 40 60 80 100 



(d) 

500 
S 300 



|T 100 - 

-100 



200 



(e) 



10 F 



CL 



1<T 



400 600 

Event Number 



800 




1000 10 2 10 ° 1 ° 2 10 * 10 ' 



FIG. 1: The difference between the activity patterns predicted by a Poisson process (top) and 
the heavy tailed distributions observed in human dynamics (bottom), (a) Succession of events 
predicted by a Poisson process, which assumes that in any moment an event takes place with 
probability q. The horizontal axis denotes time, each vertical line corresponding to an individual 
event. Note that the interevent times are comparable to each other, long delays being virtually 
absent, (b) The absence of long delays is visible on the plot showing the delay times r for 1,000 
consecutive events, the size of each vertical line corresponding to the gaps seen in (a), (c) The 
probability of finding exactly n events within a fixed time interval is V(n; q) = e~ qt (qt) n fn\, which 
predicts that for a Poisson process the inter-event time distribution follows P(t) = qe~ qT , shown 
on a log-linear plot in (c) for the events displayed in (a, b). (d) The succession of events for a heavy 
tailed distribution, (e) The waiting time r of 1,000 consecutive events, where the mean event time 
was chosen to coincide with the mean event time of the Poisson process shown in (a-c). Note the 
large spikes in the plot, corresponding to very long delay times, (b) and (e) have the same vertical 
scale, allowing to compare the regularity of a Poisson process with the bursty nature of the heavy 
tailed process, (f) Delay time distribution P(t) ~ t~ 2 for the heavy tailed process shown in (d,e), 
appearing as a straight line with slope -2 on a log-log plot. The signal shown in (d-f) was generated 
using 7 = 1 in the stochastic priority list model discussed in the Supplementary Information. 
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FIG. 2: Heavy tailed activity patterns in email communications, (a) The distribution of the 
time intervals between consecutive emails sent by a single user over a three month time interval, 
indicating that P(r) ~ r^ 1 (the solid line in the log-log plot has slope -1). While the exponent 
differs slightly from user to user, it is typically centered around a = 1. (b) The distribution 
of the time taken by the user to reply to a received message. To determine t w we recorded the 
time the user received an email from a specific user, and the time it sent a response to it, the 
difference between the two providing t w . For consistency the figure shows the data for the user 
whose interevent time distribution is shown in (a). The solid line in the log-log plot has slope 
-1. (c) A scatter plot showing the waiting time (r w ) and the size for each email responded to by 
the user discussed in (a,b), indicating that the file size and the response time do not correlate, 
(d) Scatter plot showing the number of emails received and sent by 3,188 users during a three 
month interval. Each point corresponds to a different user, indicating that there are significant 
differences between the number of received and responded emails. The dashed line corresponds to 
ni n = n ou t, capturing the case when the classic&B queuing models also predict a power law waiting 
time distribution (see Supplementary Information), albeit with exponent a = 3/2. 
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FIG. 3: The waiting time distribution predicted by the investigated queuing model. The priorities 
were chosen from a uniform distribution cc, 6 [0,1], and we monitored a priority list of length 
L = 100 over T = 10 6 time steps, (a) Log-log plot of the tail of probability P(t) that a task 
spends r time on the list obtained for p = 0.99999, corresponding to the deterministic limit of the 
model. The continuous line on the log-log plot correspond to the scaling predicted by ([%]). having 
slope -1, in agreement with the numerical results and the analytical predictions. The data was 
log-binned, to reduce the uneven statistical fluctuations common in heavy tailed distributions, a 
procedure that does not alter the slope of the tail. For the full curve, including the r = 1 peak, 
see Fig. 3 in the Supplementary Information, (b) Linear-log plot of the P(t) distribution for 
p = 0.00001, corresponding to the random choice limit of the model. The fact that the curve 
follows a straight line on a linear-log plot indicates that P(t) decays exponentially. 
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